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RESEARCH MEMORANDUM 

AN EMPIRICALLY DERIVED METHOD FOR CALCULATINC PRESSURE 
DISTRIBUTIONS OVER AIRFOILS AT ' SUPERCRITICAL MACH 
NW^BERS AND MODERATE ANCLES OF ATTACK 
By Gereld E. Nltzberg and Loma E. Sluder 

SUMMARY 

On the basis of an analysis of experimental airfoil data 
a procedure is developed for calculating pressure distributions 
over airfoils at supercritical Mach numbers. A number of 
factors limit application of the method to moderate angles of 
attack. Only the forward portion of the airfoil, where 
Interaction between local bovindary-layer thickness and local- 
pressure distribution is negligible, is considered in detail. 
Use of the method developed merely requires a knowledge of 
the low-speed or potential theory, pressure distribution, 
and the shape of the airfoil. Calculated and experimentally 
measured high-speed pressure distributions are compared for 
a wide variety of airfoils at different angles of attack 
and Mach n-umbers. 
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lOTRODUCTION 


Air, when flowing past moderately thick airfoil sections 
at seven- or- eight-tenths the velocity of sound, attains local 
velocities in excess of sonic. Analysis bf such mixed subsonl 
and supersonic types of flow presents difficulties which are 
greater than those of either purely ' subsonic or supersonic 
flow. 

A general compressibility correction for pressure coef- 
flclents has been derived by Karman and Tslen (reference 1) 
which, when applied to potential theory pressure distributions 
gives the distribution over an airfoil in any desired wholly 
subsonic flow. In the regime of purely supersonic flow, 
Prandtland Meyer (reference 2) derived a theory of the 
pressure distribution over .a curved siu?face- by considering 
the effects on an initially semi-infinite uniform flow at 
sonic velocity Then if is deflected around a -corner. If 
was found that the local supersonic Mach n-umbers attained' by ' ’ 
the- stream are a function only of. the total angle through 
which the stream is turned. This theory can be used to 
obtain the pressure distributions over' airfoils at supersonic 
Mach numbers. The supersonic flow region in the vicinity of 
airfoils at high subsonic free- stream Mach numbers is limited 
in extent so that the Prandtl-Meyer theory cannot be applied 
directly. As a- result of mathematical difficulties, a 
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readily applicable solution of the problem- of supersonic 
flow of limited extent has not been obtained as -yet^ but 
the importance of the problem makes it expedient to resort 
to an empirical analysis. . 

Pressure distributions over a number of airfoils at high 
Mach numbers have recently been obtained in the Ames 
1- by 5-l/2-foot high-speed wind tunnel. The airfoils are of 
the conventional sections NACA OOI 5 , 23015 , Uif-15 I|Ji.l2, 

and the low-drag sections KACA 652-215 (a = O.5) and 66,2-215 
(a=0.6). These pressure distributions were obtained at 
Reynolds numbers of approximately 2,000,000 and are considered 
to be accurate representations of free-^alr results up to ?£ach 
numbers of O. 8 IO for moderate an'gles of attack. The effects 
of Reynolds number variation on high-speed pressure distri- 
butions are not known. Therefore, it is not possible to 
estimate what restrictions are Imposed on the generality of 
an analysis based on these pressure measurements as a result 
of the moderate test Reynolds number. 

Experimental section drag coefficients were obtained by 
wake surveys simultaneously with the pressure distributions. 

At moderate angles of attack, the drag coefficients show no 
appreciable variation with Mach number imtil the local 
velocity of sound is exceeded at some point on the airfoil 
surface. The free- stream Mach number at which local sonic 
velocity first occurs is the critical Mach number. Above 
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the critical Mach number, there is a more or less marked 
increase in the airfoil drag coefflolont. Since a drag 
Increase is Indicatiye of a thickening boundary layer, it is 
anticipated that there are considerable boundary- layer 
effects at supercritical speeds. It would thus appear 
possible to start an analysis of the experimental pressure 
distributions by separating the viscosity and compressibility 
effects. In reference it is shown that the differences 
between airfoil section charsoterl sties in potential flow and 
in real viscous fluids can be accurately treated by taking 
account of the boundary layer. These differences were found 
to be chiefly embodied in effective changes in the airfoil 
shape and angle of attack. It thus might be hoped that the 
compressible supercritical counterpart of potential theory 
could bo obtained by correcting the experimental pressure 
distributions for the boundary- layer viscosity effects. 
Reference 1;. has treated the growth and stability of the 
laminar boiuidary layer at largo Mach n^ambers; but a survey 
of the literature reveals that boundary- layer changes through 
and behind shock waves have not yet been successfully analyzed. 
Moreover, there are a number of indications of a strong inter- 
action between the shock wave and boundary layer which 
eliminates the possibility of separating viscosity and 
compressibility effects in the vicinity of the shook wave. 

In spite of. this limitation, it still seems feasible to 
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separate gross viscosity effects, such as effective anglc- 
of-attack changes, from the compressibility effects. 

The initial step in the procedure used in the present 
report was the determination of all the known airfoil 
parameters, that is, the directly determinable properties 
of the airfoil and the flow field. These parameters were 
correlated with that property of the pressure distribution 
under consideration to find the ones which seem relevant. 

The result of such a procedure is not rigorously correct 
because it is based only on combinations of those parameters 
which are available. In other words, if some other parameter 
could be obtained, a combination of parameters might be found 
which would give equal or better correlation with the experi- 
mental measurements. Thus, the acceptance of the results of 
this report is dependent upon the rationality of the analysis 
based upon those parameters which are taken to be of prime 
importance . 


ANALYSIS 


At supercritical Mach numbers, where the flow in the 
vicinity of an airfoil is both subsonic and supersonic, it 
is to be expected that the relative Importance of local 
airfoil-shape parameters will be different for these two 
types of flow. Therefore, these two regimes will be 
considered separately. 
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Porwerd Subsonic Region 

A procedure for handling the forward subsonic region of 
the pressure distribution is suggested by analogy with the low- 
speed theory of viscosity effects on airfoil, characteristics. 
This theory shows that the pr.osence of a thick boundary layer 
results in an effective change in the angle of attack and 
camber- line.' The relatively large drag coefficients at . 
supercritical Mach n-umbers are indicative of thick boundary 
layers. Therefore, it is to be; expected that there is an 
effective anglo-of-attack change at supercritical Mach numbers. 

A study of experimental- data at supercritical Mach 
numbers reveals that, for the airfoils tested, which were of 
12- and 15-percont thickness, the subsonic .region is of very 
limited chordwise extent on at least one surface at most 
angles of attack, ' In fact, at large Mach niimbcrs the subsonic 
portion of the pressure distribution over both airfoil surfaces 
is greater than a few percent of a chord length at only 
moderate angles of attack, that is, angles from 0° to 1].°, 

So in order to handle the subsonic regime, the' analysis will 
have to be limited to such moderate angles of attack. The 
primary Importance of studying the subsonic portion of the 
pressure dlstrlbutiori is as a means for predicting the chord- 
wise location of the point at which sonic velocity is attained. 
In order to obtain a precedure which can -bo applied to an 
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arbitrarj airfoil section, it is necessary to relate super- 
critical pressure distributions to low-speed pressure 
distributions. To do this, it would seem possible to 
apply the Karman-Tsien theory to low-speed pressure 
distributions at angles of attack corrected for the high- 
speed boundary- layer effects. However, the experimental 
observations which are to be discussed directly indicate 
that in the- itnnediate vicinity of the airfoil stagne-tion 
point the variation of pressure coefficient with Mach number 
does not seen to follow the Kam^n-Tsien theory. An 
examination of subcrltical pressure distributions, at angles 
of attack for which the drag coefficient remained constant 
with the increasing Mach number, revealed that at 2|-- and 
5-percent- chord stations the experimental pressure- coeffi- 
cient variation with Mach number was less than that 
predicted by the Karman-Tsien theory, while the variation 
aft of the. 10-percent-chord station v:as in good agreement 
v/lth the theory. At supercritical Mach numbers there are 
other factors besides compressibility which affect the 
variation with Mach number of pressixre coefficients in the 
forward subsonic portion of the airfoil so that it was not 
feasible to investigate the observed deviation from theory 
in detail. Lacking a more appropriate compressibility 
correction, the present analysis will apply the Karman-Tsien 
theory even near the airfoil leading edge. 
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In low-speed viscous theory, an effective angle-of-attack 
change is considered to occur when the presence of a boundary 
layer alters the flow past an airfoil in a manner which is 
similar to that which would arise from a shift of the airfoil 
trailing edge from its actual physical position to an 
effective position at the center of the total trailing- edge 
boundary layer. Thus, to find the effective angle-of-attack 
change, it is only necessary to know the trailing- edge 
boundary-layer thickness for each surface of the airfoil. At 
low speeds the drag coefficient is a known function of boundary- 
layer thickness and pressure coefficient at the trailing edge. 
However, this functional relation is not known for super- 
critical speeds. Moreover, at supercritical Mach numbers, 
drag has two components - losses through the shock wave and 
within the boundary layer. It is not to be expected that tho 
shock-wave wake would operate like the boundary layer, in 
changing the effective angle of attack. Therefore, the 
traillng-edge boundary- layer thickness cannot bo directly 
determined at supercritical Mach numbers from the known 
pressure distribution and drag coefficient. 

That portion of the experimental pressure distributions 
ahead of the sonic point was reduced to low speeds by using 
the Karman-Tslen theory, VJhen these adjusted pressure 
distributions were compared with loiv-speed pressure distri- 
butions, it was found that the upper and lower surfaces of an 
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airfoil appeared to be at different angles of attack. The 
adjusted pressure distributions over single surfaces did 
not conform exactly to the low-speed distribution for any 
angle of attackj however, over the forward subsonic portion 
of each pressure distribution, the maximum scatter 
corresponded .to a spread of a few tenths of a degree. This 
sim.ilarity indicates that, as a first approximation, thc- 
total boundary- layer effect can be treated as an angle-of- 
attack chan,ge. A precise evaluation of boundary- layer 
effects would require separating the angle-of -attack 
effect and the shape change. There is no available method 
for calculating boundary-*layer changes in the vicinity of a 
shock wave so the present analysis is limited to the 
previously mentioned first approximation. 

The effective angle of attack for each airfoil surface 
was calculated from the experimental pressure distributions 
for which only one surface had a region of. supersonic- flov/. 

For the moderate angles of attack considered, the supersonic 
region always occurred on the upper surface of the -airfoil. 

The values of the effective angle-of-attack change, occasioned 
by the presence of a shock on the upper surface Aa(^»S,S.) 
were plotted for each surface as a function of the Increment 
of the free- stream Mach number above the critical Mach 
number M-ItQj.. Values were plotted for the upper surface in 
figure 1 and those for the lower surface in figure 2, In these 
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figures it is seen that for the wide variety of airfoil shapes 
and differences in free- stream Mach number and geometric angles 
of attack, the maximum dispersion of the points for any value 
of M - Is only about lo. Therefore, a curve faired through 
the mean values of the experimental points can bo used to 
predict the effective angle-of-attack change for a given value 
of M - to within about l/2°. . Although this' approximation 

is adequate for predicting the sonic-point location, an 
attempt was made to determine whether this . dispersion was a 
function of some other parameter, A parameter which seemed to 
be of considerable importance was dlJ^jp/da , It v/as noticed 
that the majority of the cases for which Aa(U.S,S.) values 
were low corresponded to configurq^tlons having values of 
dMcr/da greater than 0,025. Further evidence of the 
importance of this parameter is presented in figure 5 whore _ 
the Increment of drag coefficient above the low- speed drag 
coefficient Op - is plotted as a fxinction of M - Mqj,, 

for the configurations considered in figures 1 and 2, Figure 5 
shows that, for all but one conf Igioration, a curve can be 
drawn separating the points having large and small values 
of dMcj./da . A method, . based on this parametei) has been 
developed for predicting the effective angle of attack with 
greater accuracy; however, the theoretical critical Mach 
numbers are not in sufficiently close agreement, with 
experimentally measured values to justify the use of this 
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more involved method. 

The critical Mach number of an airfoil at a given geo- 
metric angle of attack, can be found by applying the Ka'ra^n- 
Tsien compressibility correction to the airfoil peak-pressure 
coefficient obtained from potential theory. This procedure 
has been used in reference 5 * It is found that the variation 
of critical Ilach number of each surface of an airfoil with 
angle of attack is different for airfoils of different shapes, 
and,- even for a glvon airfoil, there are tv/o markedly different 
anglc-of- attack regimes for each surface in which this differ- 
ence in variation is apparent. One regime, usually positive 
angles of attack for lower surfaces and negative angles for 
upper surfaces, includes angles for which the peak-pressure 
coefficient is determined chiefly by the airfoil- thickness 
distribution and camber, and consequently there is little 
change of critical Mach number with angle of attack. In the 
other regime, the maxlmxim-pressure peak occurs near the air- 
foil nose as a result of the additional lift distribution, 
that is, the increased circulation resulting from pitching 
the airfoil. There is considerable uncertainty as to the 
accuracy of the calculations of critical Mach numbers at 
these angles because of the previously discussed deviation 
of the flow near the stagnation point from the Karmefn-Tsien 
theory. In this regime, uncertainty also exists for critical 
Mach numbers determined from high-speed experimental presaiare 
distributions unless a number of pressure orifices have been 
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pi". cod . along the forward 5 percent of the airfoil chord. This 
lin'.itationa places serious restriction on a significant source of 
information regarding the high-speed characteristics of airfoils. 

The experimental and -theoretical critical Mach number 
curves for* both surfaces' of the airfoils considered are shown 
as a function of angle of attack- in figure i;. As will be 
discussed directly, in the present method of analysis, a 
correction to the theoretical curves is required, 

WTien only one surface is at a supercritical Mach number 
there is an effective angle-of-attack change of the other 
surface which results in a change in the pressure distribution 
and thus in the critical Mach number of the latter surface. 
Therefore, if the critical Mach number curves are obtained 
theoretically, it is necessary to adjust the calculated 
values for the surface with the higher critical Mach number 
to account for the effective angle-of-,attack change due to 
the shock wave on the other surface of the airfoil, Cvirves 
thus derived are in reasonable agreement with experimental 
measurements . 

As was indicated previously, when the free- stream Mach 
n^umber M is above the critical Mach number of only one 
surface of an airfoil, the following procedure can be 
used to calculate the subsonic portion of the pressure 
distribution over the airfoil. The critical Mach number 
Mqp, for the moderate geometrical angle of attack under 
consideration, is computed from the low-speed pressure 
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distribution by means of the Kdrm^n-Tsien theory. Curves of 
the Kdrmdn-Tsien compressibility correction are given in 
figure 5* Then having the value of li - the curve of 

figure 1 ia used to estimate the effective angle of attack 
for the surface of the airfoil over which the local region 
of supersonic flow occurs. The subsonic portion of the 
pressure distribution over this surface of the airfoil at 
the f ree-stream i^sch number under consideration is obtained 
by applying the Kariran-Tsien compressibility correction to 
the low-speed pressure distribution corresponding to the 
effective angle of attack. The flow over the other airfoil 
surface is entirely subsonic so that the pressure distri- 
bution can be estimated directly. The effective angle of 
attack of this surface is obtained from the curve of figure 2 
by using the same value of M - Hcp as used in figure 1. 

The previously mentioned procedure can be generalized to 
apply to the ca<=e of airfoils at Mach nximbers above the 
critical speed of both surfaces if the assumption is made 
that the effects of the supersonic region on each surface 
are additive. In this case, it is necessary to obtain the 
critical Mach number for each surface of the airfoil. Then 
using the upper-surface Mach number increment above its 
critical (M - Mcr^yA effective angle-of-attack changes, 

due to the upper- B^ar face supersonic region of flow, for the 
upper sinrface AaCU.S.S.)^ and lower surface Aa(U.S.S.)L 
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of the airfoil are obtained from flares 1 and 2, respectively. 
Similarly, for the lower-surfao© Increment above Its critical 
(M - ^^crjj) effective arxgle-of-attack changes, due to the _ 

lower- surface supersonic region of flow, for the lower surface 
Aa(L.S.S.)j^ and upper surface Ao(L.S.S.)y are obtained from 
figures 1 and 2 , respectively. The total effective angle-of- 
attack change Is then Aa(U.S,S, )y - Aa(L.$.vS.)y for the upper 
surface, and Acc(U,S.S. - Aa{L.S,S.)L for the lower surface. 
The limited experlii^H^iftl to derive the afore- 
mentioned general proceclui*® it of dubious accuracy at 

Mach numbers above M - of 0 , 20 , Below this limit the 

maximum difference betweeh the calculated angle of attack of 
a surface, predicted by this procedure, and that obtained from 
the experimental pressure-, distribution was about l/2°. 

The preceding method can be used to predict the chordwlso 
location of the point at which sonic velocity is first attained. 
The pressure coefficient corresponding to sonic velocity at a 
given free- stream Mach number may be obtained from figure 5 * 

The movement of the sonic point with Mach number has been 
predicted using this procedure for the airfoils for which 
experimental pressure distributions are available. The method 
is applicable for the normal operating range of high-speed 
airplanes, that is, angles of attack from 0 ° to I4-® and Mach 
n\imber increments less than 0,20 above the critical, Jn 
this range the difference between the calculated value and 
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the experimentally measured sonic point v/as less than a few 
percent of a chord length. 

Region of Increasing Supersonic Local Mach Numbers 

The Earman-Tsien compressibility correction is strictly 
limited to subsonic regions of the pressure distribution. 

The only available procedure for supersonic flows is that of 
Prandtl and Meyer, which presents the local Mach number as a 
function of the difference between the local surface slope 
and the slope at the sonic point. The local Mach number is 
uniquely related to the pressure coefficient for a given 
free-stream Mach number. This relationship is shown in 
figure 6 . To study the supersonic portion of. the airfoil 
pressure distribution at moderate angles of attack, local 
Mach number was plotted as a function of the change in slope 
of the surface from that at the sonic point. It was found 
that the variation was similar to the Prandtl-Meyer theory 
but different configurations achieved varying percentages 
of the theoretical local Mach number rise above unity 
corresponding to a given slope change. The Prandtl-Meyer 
theoretical increase in local Mach number for expansive 
deflection of the stream from its direction at sonic 
velocity and several percentages of this theoretical 
variation are presented in figure ?, For a given airfoil, 
angle of attack, and Mach number, the local Mach numbers 
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attained conformed closely to a unique percentage of the 
Prandtl-Meyor theoretical results. This failure to obtain the 
full theoretical rise of local Hach number seems reasonable 
because the Prandtl-Meyer theory Is baaed on the assumption 
of a seml-lnflnl te ■uniform sonic flow deflected around a 
curved surface; whereas, airfoils at supercritical Mach 
numbers have a supersonic region of very limited extent. The 
extent and uniformity of this supersonic region depends on 
such factors as the free- str earn l^ach nvimber, the critical 
Mach n\Amber, and the chordwlse location of the sonic point. 

One variable ■ chosen was “Percent of Prandtl-: Moyer” and, 
this was correlated with- all known airfoil parameters which 
seemed relevant. These parameters included H - MQj,,-the 
surface curvature at the sonic point and also at the point at 
which a local Mach n'umber of 1,2 occurs, the slope of the 
airfoil surface at the sonic point re3^tlve to the free- 
stream direction and also relative to the airfoil chord line, 
the chordwlse location of sonic velocity and local Mach number 
of 1*2, and others. It was found that the two parameters 
which gave the best correlation with percent of Prandtl-Meyer 
were H - M^p and the slope of the surface at the sonic point 
relative to the airfoil chord line. Measuring slope relative 
to the airfoil chord line seems strange at first consideration 
but, since for the configurations considered, the stagnation 
point very nearly coincides with the forward end of the chord 
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line, there seems to he some physical basis for using the 
airfoil chord as a reference lino. The data which were 
analyzed Included ca.ses for w'hich the Included angle 0 3 
between a tangent to the surface at the sonic point and the 
airfoil chord line varied from 5° 50°. yJlthin this range 

it was observed that for equal values of M - the percent 

of Prandtl-Meyer attained decreased with increasing 83^ It 
should be mentioned that increasing values of 63 correspond 
to movement of the sonic point toward the airfoil leading 
edge. The other variation noted was that increasing M - 
for constant values of 83 resulted in increasing the 
percent of Prandtl-’leyer which was attained. Curves based 
on the experimental data are presented in' figure 8, The 
maximum difference between experimental values and these 
curves corresponded to a 5~P®^cent shift in percent of 
Prandtl-Meyer, 

Once the slope distribution over the airfoil and the 
sonic-point location are determined, the supersonic portion 
of the pressure distribution can be calculated as follows: 

The percentage of the Prandtl-Meyer theoretical curve which 
should be used is determined from figure 8, The local Mach 
numbers reached at any chordwise station are obtained from 
the particular percentage curves shown in flgTire 7* A'hen 
the curves of figure 6 are used, these local Mach numbers 
can be converted directly to pressure coefficients S. 
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Of course, this analysis is based on the use of perpent. 
of Prandtl-TIeyer as a parameter and thus can only be considered 
to be empirically derived. 

Termination of the Supersonic Regime 

There is nothing in the foregoing analysis to give any 
clue as to the chordwise point of termination of the super- 
sonic region. It was believed that there was a very direct 
relation between the stability of the boundary layer and the 
location of the shock wave. To check this hypothesis, 
pressure distributions over a carefully cleaned NACA [jl|.12 air- 
foil were obtained at a moderate angle and various Mach nvimbers 
in the Ames 1- by 3“l/2-foot high-speed wind tunnel. Then 
grains of carborundum were attached between the 5“ 10- 

percent chordwise stations on the airfoil, and pressure 
distributions at the same angle and Mach numbers were, obtained. 
The pressixre recovery, corresponding to the shock wave, was 
much more abrupt but only slightly further forward for the 
roughened airfoil than for the smooth airfoil. The pressure 
distribution in the region of increasing supersonic local 
Mach numbers was the same for the two airfoil surface conditions. 
It is not feasible to make a quantitative analysis of this 
phenomenon since there are no available theories of laminar 
separation, transition from laminar to turbulent boundary- 
layer flow, or turbulent separation, in supersonic flows. 
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Even more Important, there are no quantitative results for 
the interaction between such boundary- layer changes and shock 
waves. For the pressure distributions available on airfoils 
of 12- and 15-percent tliickness, for moderate angles of 
attack at Mach numbers between O.O5 and 0.20 above the 
critical, it was noticed that the point at which the local 
supersonic Mach numbers ceased to increase corresponded 
rougjhly to that point, on the surface under consideration, 
at which the surface was tangent to the free- stream direction 
It is not immediately apparent why this geometrically deter- 
mined parameter should so Influence the flow past the airfoil 
There is usually ho precise location of a shock wave 
discernible from the pressure distributions. -Instead of an 
abrupt discontinuity there is a gradual pressure recovery. 
These pressure distributions were obtained by means of a 
liquid manometer so that they may represent a time average 
pressure ■distribution of an unsteady flow condition. 
Alternatively, the gradual pressure recovery observed may be 
a result of a thick boundary layer. This bo-undary layer may 
soften the shock -as discussed in-reference 6-. Also,' the 
abrupt pressure change might be masked by the boundary- 
layer thickness, but this seems ixnllkely since, as has 
been previousl'j mentioned, an abrupt pressure recovery was 
observed on the roughened NACA i|l|.12 airfoil. It is expected 
that, since airfoil surface roughness causes marked changes 
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in the nature of the pressure recovery which accompanies 
transition from supersonic to subsonic flov/, Reynolds n-umber 
might have important effects on the shock wave. In all tho - 
configurations analyzed it was found that the position .pn the 
airfoil surface at which the local Mach number had dropped to 
unity could be estlmatedby applying the Karman-Tsien 
compressibility correction to the potential theory pressure 
distribution. It was the-refore decided to terminate the 
supersonic portion of tb© pres sure- distribution calculations 
at the chordwlso station for which a local Mach number of 
unity is predicted. 

The experimental pressure dlstributlcns indicate a more or 
less gradual pressin’o recovery over the region from tho chord- 
wise station at which the surface ivS tangent to the free- 
stream direction to that station at which a local Mach number 
of unity occurs. Valien tho local pressxare coefficients for 
these two. stations are connected by a straight lino, it is 
noticed that the local portion of the “experimental pressure 
distribution lies near this line or between this lino- and the 
calculated super.sonic pressure distribution. The presence of 
this region of decreasing supersonic local Mach numbers is 
believed to be due to Interaction between the local boundary 
layer and the shock wave. 
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Estimation of Pressure Distribution Aft 
of Supersonic Region 

Over the rear portion of an airfoil, behind a shock wave, 
the bo\mdary layer is probably quite thick. It is to be 
expected that a thick boundary layer could have considerable 
effect on the local pressure distribution. However, there 
are no available methods for calculating the boundary- layer 
thickness behind a shock wave. Therefore, in the present 
report the subsonic portion of the pressure distribution 
over the rear of the airfoil will be approximated by 
applying the Karman-Tsien compressibility correction to the 
local potential theory pressure distribution. 

DISCUSSION 

The primary concern of the preceding analysis is the 
prediction of supercritical pressure distributions over the 
forward portion of airfoils. The detailed analysis is 
restricted to this portion of the airfoil where local 
boundary- layer thickness does not markedly alter the local 
pressure distribution. Before comparing calculated and 
experimental pressure distributions, some of the general 
problems which arise v/ill be reviewed. 

In comparing the low-speed or potential theory pressure 
distribution to those at moderately high Mach numbers, it is 
found, as would be expected, that the Karman-Tsien theory 
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is Inaccurate near the nose of the airfoil. There is no 
available procedure for treating this region. Therefore, in 
the preceding analysis the assumption is made that the TC.rrmah- 
Tsien theory provides an adequate approximation for pressure 
coefficients at points close to the airfoil leading edge. 

However, the errors introduced by the assumption are of secondary 
Importance at moderate angles of attack. An investigation is 
needed to determine how pressure coefficients change with Mach 
number at the nose of the airfoil. 

The airfoil- surface slope changes rapidly in the vicinity; 
of the airfoil nose. When the sonic point occurs in this region, 
there may be considerable error in the determination of the 
surface slope at the sonic point. Fortunately, the present 
method compensates for errors in the sonic-point slope. For 
example, if too large a sonic-point slqpe is read for a given 
value of H - M^p, too small a percentage of Prandtl-Meyer 
theory would be obtained from figure 8, 'linco local Mach number 
attained varies with both the surface- slope change from the sonic 
point and the percent of Prandtl-Meyer, and an excessive value 
of one results in a diminished value of the other, it is seen 
that these two effects tend to counterbalance each other, except 
in the immediate vicinity of the sonic point. 

The present method can be expected to give rather good 
accuracy in determining the peak-pressure coefficients and the 
maximum load obtained on the forward portion of the airfoil, 
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Hov/ever, a factor which may contribute to differences between 
calculations and experimental measurements for airplane wings 
is the fact that the supersonic portion of the pressure distri- 
bution depends upon the airfoil- surface slope. As a result of 
practical construction difficulties, there are discrepancies 
in the actual wing profile v/hich give rise to uncertainties 
regarding the surface slope. 

A.ppll cation of Method 

In order to evaluate the utility and limitations of the 
procedure which has been developed, this method has been used 
to calculate the various configurations for which experimental 
results are available. As has been pointed out, the develop- 
ment of the method limits its application to moderate angles 
of attack and to Mach numbers less than 0.20 above the critical 
Mach number. For angles of 0°, 29, and two supercritical 

Mach numbers have been considered for each airfoil. The calcu- 
lated characteristics of all the configurations are tabulsted 
in table I. 

For incompressible potential theory, it is usual to obtain 

the velocity distribution rather than the pressure distribution. 

The additional lift due to pitching the airfoil is directly 

additive to the basic velocity distribution so that the velocity 

distribution is particularly convenient for treating cases for 

which several angles of attack are considered. Since velocity 

distribution is so easily modified, it was decided to use the 
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pressure coefficient which, at low speeds, is the square of 
the ratio of the local velocity and the free- stream velocity. 
In figure 5 the Karman-Tsien theoretical variation of pressure 
coefficient S with Mach member is presented. The supersonic 
theory of Prandtl-Heyer is in terms of the local Mach mmiber. 
Figure 6 gives the pressure coefficient 5. ..corresponding to_..a 
local Ilach number as a function of free- stream Hach number. 

As has been previously mentioned, it was noticed that the 
termination of the supersonic regime of rising pressure coef- 
ficients occurs approximately at the point where the airfoil 
surface is parallel to the effective free-stream direction, In 
the calculated pressure distributions this point is connected 
with the rearward sonic point by a straight lino. 

As an example of pressure-distribution variation with 
Mach number for a fixed geometric angl^ of attack the NACA 
652-215 airfoil at 2 ° is considered. Figure 9 shows a com- 
parison of the calculated pressure-distribution curves and 
the experimental points at Mach numbers of O.'JOO, O.725, 

0.750, 0,775, 0.800, and 0.825. The present method has a 
lower limit of applicability at - Her of O.O5 because of 
the very limited supersonic region for such slightly super- 
critical Mach numbers. An approximate value of the pressure 
distribution in this regime is obtained by applying the 
Ka'rman-Tsien compressibility correction for the entire 
pressure distribution. At Mach numbers of 0,775 0,8 00 
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the lower surface is only slightly supercritical and applying 
the Ka'msCn-Tsien theory gives inaccurate values for the peak- 
pressure coefficients. Both surfaces are well above their 
respective critical Mach numbf-rs at O .825 so the present 
method can bo applied to both surfaces, The large drag 
coefficient at this T^ach number Indicates a thick boundary 
layer at the rear of the airfoil which may be the cause of 
poor accuracy of Ksfrmah-Tsien theory in this region. Through- 
out the Mach number range the calculations give a good 
prediction of the peak-pressure coefficients and of the load 
over the forward I 4 .O percent of the chord length. 

Now consider the same airfoil at different angles of 
attack. At and the marked difference in the type of 
loading for the low and high supercritical Mach nu’^bers is 
correctly represented by the calculations. 

In figure 10, it is seen that the method works well for 
the NACA 66,2-215 airfoil at 0° and 20, As a result of the 
small leadiag-edge radius of this airfoil the upper-surface 
critical speed at I|.o is determined by a pressure peak near 
the nosej therefore, the critical speed cannot bo determined 
accurately,. However, this is a borderline case in that the 
forward pressure peak is not very sharp, so the critical Mach 
number v/as based on the farthest forward experimental pressure 
point. The calculated pressure distributions based on this 
critical Mach number are presented to show the magnitude of 
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the error that can be expected for such borderline cases. 

■ For the NACA i(i4.12 airfoil at 0°, 2°, and i|.° the change 
in typo of loading with Mach number is accurately predicted. 

There was some inaccuracy in the construction of the 
NACA. 0015 airfoil as can be seen in figure 12 from the 
difference between the experimental pressure distribution for 
the two surfaces at 0°. At 2P for » 0.725 and at [{.o for 
both Mach numbers shown, the upper-surface shock is very far 
forward. The start of the marked pressure recovery is still 
at the point at which the airfoil surface is tangent to the 
froe-stream direction. 

The comparisons of experiment and theory for the NACA 
23015 are presented in figure 15. At'Q'^, 2“^, and ii.*^ tho shock 

occurs far forward, but the load distribution and pressure 
peaks are well represented. 

It is interesting to notice the similarity between the 
results for the NACA I4.l4.i5 airfoil in figu:i:*e li|- and the NACA 
ljll.12 in figure 11, The values of M * used in figure 1J4. 
were chosen to be as close as possible to those for the 
corresponding angles of. attack in figure 11, The agreement 
between calculation and experiment in the vicinity of the 
nose is not as good for the NACA i|J4.15 airfoil as for the 
NACA I4I4.I2 airfoil because of the difficulty arising as a 
result of the large leading-edge radius of the former airfoil. 
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CONCLUDING REMARKS 


The procedure which has been developed provides a satis- 
factory prediction of the supersonic portion of the super- 
critical pressure distributions for airfoils at moderate 

* 

angles of attack. The parameters on which the procedure 
is based were limited to those which were determinable 
from potential theory pressure distributions and the 
geometric properties of the . airfoils . This method, therefore, 
cannot bo considered definitive. It is primarily a useful 
means for predicting maximum loads and peak-pressure coeffi- 
cients, Further experimental investigation of conditions in 
the vicinity of the airfoil nose at high Mach numbers should 
permit an extension of the method to a wider angle-of-attack 
range , 
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